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FOKEWOBD 


The motivation for the research reported -in this document was 
to delineate the adverse effects of wind shear on the landing 
flight phase of aeronautical systems. Once these effects are 
known, relative to the total wind environment, it is possible to 
establish operational wind shear requirements and limits for observ- 
ing and reporting low level wind shear. To assess the effects of 
wind shear, or to grade wind shears relative to their effect on 
aeronautical systems, a criteria must be established in the context 
of aeronautical sysfem performance parameters. In- this study, the 
degree to which a given shear environment adversely effects the 
landing flight phase of aeronautical systems was assessed in 
terms of the departure of the landing touchdown point from the 
touchdown point that would have occurred in the absence of wind 
shear. A variety of aircraft types and a broad selection of 
wind environments were selected for the analysis. The selected 
wind environments encompass a significant number of low level 
wind situations that would be encotmtered during the life of an 
operational aeronautical system. The effects of the wind environ- 
ments on the aircraft touchdown point are presented in terms of 
properties of the selected flow fields. A number of new conclu- 
sions resulted from the study relative to how the details in 
the wind profile can effect the landing flight phase. It is 
believed that these results can have significant Implications 
relative to the aeronautical safety aspects of the landing 
problem. 

This research was conducted by the University of Dayton Research 
Institute for the National Aeronautics and Space Administration, 

George C. Marshall Space Flight Center, Huntsville, Alabama, 

under the technical direction of Dr. George H. Fichtl and Dr. Stephen W. 

Winder of the Aero-Astrodynamics Laboratory. The support for this research 

was provided by Mr. John Enders of the Aeronautical Operating 
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NASA Headquarters. 
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INTRODUCTION 


Wind shear is an important consideration in the landing of aircraft 
and aerospace vehicles. As an aircraft descends its glide slope, a sudden 
change in the h orizontal wind will instantaneously effect the velocity of the 
aircraft relative to the air mass.* If the shear is such that the relative 
velocity of the aircraft increases, the lift force will increase and the air- 
craft will tend to rise above the glide slope. If the shear causes a sudden 
decrease in the relative velocity, the aircraft will respond by falling below 
the glide slope and a potentially hazardous condition could result. 

Several reports have been published which link short and long 
touchdown to a sudden wind shear occurrence during final approach (Ref- 
erences I and 2). Recent accident reports have found wind shear to be at 
least a contributing cause to several accidents (Reference 2). In addition, 
it is believed that wind shear has been responsible for many other accidents 
though it remained undetected at the time (Reference 3). 

The problem of quantitatively defining the effect of shear of given 
magnitude on an aircraft during descent has not been completely resolved. 
Noteworthy studies that have investigated wind shear and/or turbulence 
during landing include References 4, 5, 6, 7, 8, and 9. The study under- 
taken by the University of Dayton Research Institute (UDRI) was designed 
to provide answers to three specific questions: a) what shape of wind 
shear profiles are most critical to aircraft landing, b) which type of air- 
craft are most responsive to shear, and c) what meteorological parameters 
relate to those wind shears that provide critical landing problems. The 
UDRI provided answers to these questions by a digital simulation model 
for an aircraft landing in various wind profiles. The simulation model 
is used by first calculating the touchdown point for a conventional -type 


* The effect of vertical air motions is not considered in this report. 



aircraft trimmed on an initial glide slope of 2. 7 degrees (7. 0 for a STOL 
aircraft) in a constant wind field descending from an altitude of 300 feet. The 
landing simulation is then repeated for a wind shear profile input with initial 
trim conditions determined at 300 feet for the wind velicity at that altitude. 

By determining the deviation in touchdown point from the constant wind case, 
the effect of the shear on aircraft touchdown can be determined. This simu- 
lation model assumes a fixed stick model with no pilot or autopilot control. 

The wind shear profile is defined in the surface boundary layer 
according to similarity theory by the surface roughness length, Z^; the 
zeroplane displacement, d; the surface friction velocity, u*; and the sta- 
bility parameter, Z/L. The next two sections describe in further detail 
the landing simulation model and the wind profiles used in the study. 

Throughout this report, English units are used to describe the 
aircraft-related quantities while metric units are used to describe the 
meteorological quantities. On all figures, a dual system consisting of 
both sets of units is used. The use of both unit systems was necessitated 
to conform to the conventional systems used by the aircraft and meteoro- 
logical communities. 


2 


AIRCRAFT LANDING MODEL 


The aircraft trajectory model employed in this study was derived 
based on the following assumptions. 

a) The earth is flat and nonrotating. 

b) The acceleration of gravity is constant (32. 2 ft/sec^). 
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c) Air density is constant (0. 002375 slug/ft ). 

d) The airframe is a rigid body. 

e) The aircraft is constrained to motion in the vertical plane. 

f) The aircraft has a symmetry plane (the x-z plane). 

g) The mass of the aircraft is constant. 

h) Once the aircraft is trimmed, its throttle setting and elevator 
deflection angle are not changed. 

i) The aerodynamic stability derivatives are constant within the 
altitude and Mach number range experienced in this investigation. 

At the beginningofeachtrajectory (300 feet altitude, H), the aircraft is trimmed 
by determining the values of angle of attack, throttle setting, and elevator 
deflection, which will result in unaccelerated flight. The equations of motion 
are then integrated numerically by a fourth-order Runge-Kutta scheme. 

For a constant wind and no ground effects, the aircraft flies down the glide 
slope at a constant velocity until it reachs the ground. Upon introducing a 
vertically -varying horizontal wind field, the aircraft no longer adheres to 
the glide slope. The resulting deviation in touchdown point serves as a 
measure of how severe the particular wind field is to the landing aircraft. 

The influence of ground effects on the deviation in touchdown points 
between constant wind and wind shear conditions was investigated for several 
of the aircraft considered in this study. It was found that this influence was 
rather small and therefore was not included in the final analysis since 
ground effects data was not available for all of the aircraft. This is dis- 
cussed further in the section describing the analysis of the data. 
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The aircraft included in this study are the DC-8, C-135A, C-141, 
C-130E, B oeing-747, and an augmentor-wing STOL aircraft. The first three 
are representative of the medium-weight turbojet transports. The first 
two are low-wing design while the third is a high-wing design. The C-130E 
is a lighter- weight transport powered by propjet engines. The Boeing -747 is, 
of course, a large turbojet transport. The augmentor-wing STOL aircraft 
is in the same weight category as the DC-8, C-135, and C-141. The aero- 
dynamic data for this aircraft is similar to that used by NASA Ames 
Research Laboratory in their computer simulations of an augmentor-wing 
STOL aircraft. 

The equations of motion of the aircraft were derived under the 
assumptions stated above. Figure 1 shows the forces acting on the aircraft. 

Ok 

These include gravity (mf"), thrust of the engines and the aerodynamically- 

induced lift (L), and drag (D) forces. The figure shows the orientation of the 
forces with respect to the velocity vector relative to the earth (V), the velocity 
vector relative to the air mass (V-^), and the fuselage reference line (FRL) 
of the aircraft. The X axis in Figure 1 is parallel to the surface of the 
earth; the Z is perpendicular to the surface of the earth (positive downward). 

Two of the equations of motion can be derived by summing the 
forces parallel and perpendicular to V (the velocity vector relative to the 
earth) and applying Newton's Laws of Motion. The result is 

^T 

V = - g sin V + cos (6_ + a ) 

cos6 (1) 

C sin6 . 

Xj 


m 

si 

m 


4 




Figure 1. Relationship Between the Various 

Forces Acting on an Aircraft. ' 

i 
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- Y sin(5^+a) 

+ C cos6 , 

mV L 


C sin6 
mV D 


where the dot refers to the derivative with respect to time and 

g is the magnitude of the acceleration of gravity, 

V is the magnitude of the velocity relative to the earth, 

Y is the angle between V and the X-axis (the flight path angle), 
is the magnitude of the thrust vector, 

m is the aircraft mass, 

6^ is the angle between the thrust vector and the fuselage 
reference line (FRL), 

a is the angle between V and the FRL, 

q is the dynamic pressure which is equal to one half the air 
density (p) times the square of magnitude of the velocity 
relative to the air mass (V^^), q = 1/2 pV^^ 

5 is the aircraft wing area, 

6 is the angle between V and V, 

Cjj is the drag coefficient, and 

C. is the lift coefficient. 


The aerodynamic forces and the thrust from the engines exert a 
pitching moment on the aircraft. The equation describing the rotational 


acceleration is 


T T 


is the time derivative of the pitching rate (q), 
is the effective moment arm of the thrust vector, 
is the Mean Aerodynamic Chord, 
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is the moment of inertia about the symmetry plane of the 
aircraft, and 

is the pitching moment coefficient. 



Equations (1), (2), and (3) form the core of the Aircraft Landing 
Program. These three equations along with 

X = V cosy 
Z = -V sin-y 


are the six equations of motion of an aircraft constrained to fly in the 
vertical plane. 0 is the time derivative of the angle between the FRL and 
the X-axis. 


In order to evaluate the above equations at eachtime step, several 
aiixiliary equations are needed. They are 

Va = [(X-W^)2 + Z^J 
q =l/2pV^2 
-1 

6 = sin ( — ) 

a’ = 0 --y - 6 9 


a 


a ^ 


= 0 . Y 

= cosW+M 




Cj. 'C l**'" *E-\- > 
'e- Va--)- "'). 
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where "W^ is the horizontal wind speed, a* (the angle of attack) is the angle 
between and the FRL, and 6g is the elevator deflection angle. As is 
indicated above, the aerodynamic coefficients are functions of a number of 
variables. The expressions for C_ , C_, and C , are not the same for all 
the aircraft considered. These expressions along with the stability deriva- 
tive data are presented in Appendix A. 

The above set of equations comprise the aircraft model. The initial 
conditions and aircraft physical data for each of the flights simulated in 
this study are presented in Table 1. The glide slope for the conventional 
aircraft flights was 2. 7 degrees while its value for the STOL aircraft 
flights was 7 degrees. As noted in the table, the Boeing -747 was investigated 
at two landing weights. Flights of this aircraft were simulated with center- 
of-gravity locations of 15%, 25%, and 33% Mean Aerodynamic Chord. Each 
aircraft was flown in a wide variety of horizontal wind shear conditions. 

These are discussed in the following section. 
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WIND SHEAR MODEL 


The wind shear model used in the aircraft landing simulation is 
completely described in the document, "A Model of Wind Shear and Tur- 
bulence in the Surface Boundary Layer" (Reference 10). Only a cursory 
description of the model is presented here. 


The wind shear in the surface boundary layer is considered to be 
a function of surface conditions, stability conditions, and altitude. According 
to similarity theory, the mean wind speed for three of the four stability 
classifications is defined as a function of altitude by 

u = [Ln (^) + (4) 

k Zq L 

where is the surface roughness length, 
u* is the surface friction velocity, 
k is the Von Karman's constant *0.4, 

Z is the altitude above the reference level, and 
L is the Monin-Obukov stability length. 


In the unstable classification, the function ilr is given by 

Z / L 

= [ Z/L I 1- (1-18Z/L I d(-^) 

^ -'Z /L 


The stability parameter, Z/L, is related to the gradient Richardson 
number for the unstable classification by Businger's Hypothesis: 


Ri = Z/L for Ri < 0 . 

For the neutral classification, Z/L = 0 anditt(0)=0 so that the wind speed, 
given by Equation (4), is a logarithmic function of altitude. 
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For the stable condition, that is 0 < Ri ^0. 2, the function is described 
by 

♦ (Z/L)=5. 2Z/L 

The relationship between Richardsons number and Z/L is stable air 
is 

Z/L = , for0<Ri<0. 2 

Figures 2» 3, and 4 show plots of typical wind profiles for the unstable, 
neutral, and stable categories. 

For the very stable condition, when Ri^O. 2, the wind speed cannot 
be represented by Equation (4). In fact, no analytic function has been foxmd to 
adequately represent the very stable wind profiles. Under strong inversion 
conditions, turbulence tends to diminish so that the layers of the atmosphere 
become disconnected. Figures 5 and 6 show two types of profiles that are 
likely to exist under inversions. Figure 5 shows a calm below the interface 
level with a constant wind above this level. In Figure 6, a logarithmic wind 
profile is shown below the interface with a constant wind magnitude above 
this level. Several values of the interface level and the wing magnitude above 
the interface were considered in this study. 

The wind direction in the surface boundary layer can be considered 
constant with altitude except in the very stable condition. In very stable 
air, the wind direction often changes by 45 degrees or more between the 
surface and 300 feet. 

The wind input to the aircraft landing simulation program consists 
of defining the wind profile from 300 feet altitude to the surface. The 
wind magnitude has been defined for the unstable, neutral, and stable 
conditions by Equation (4). The parameters u*, Z^, and L were varied 
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Figure 5. Very Stable Wind Profiles. 



Figure 6. Very Stable Wind Profiles, 
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so as to include all reasonable wind profiles. For the very stable condition, 
wind profiles of the type shown in Figures 5 and 6 have been used. The 
wind direction is considered constant with altitude for all stability conditions 
In the landing simulation program, only headwinds and tailwinds have been 
considered with the emphasis being on the more conventional headwind 
landing case. 
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ANALYSIS OF CONVENTIONAL AIRCRAFT LANDINGS 


The initialization conditions for the simulated landings of conven- 
tional aircraft were a 2. 7-degree glide slope with the descent beginning 
at an altitude of 300 feet. (An augmented -wing STOL aircraft is discussed 
in a later section. ) This corresponds to a touchdown point 6361 feet down- 
range from where the descent begins. The aircraft is trimmed at 300 feet 
to maintain the 2. 7 -degree glide slope for the wind speed existing at that point. 
Any variation in wind speed will cause the aircratft to deviate from the glide 
slope. The deviation in touchdown point is defined as the distance between 
the actual touchdown point and the 2. 7-degree glide slope touchdown point. 

That is, if the aircraft lands at a distance X downrange from its initial 
descent point (300 feet altitude), then the deviation in touchdown point, AT, 
is 

AT = X - 6361. 

Note that a positive AT indicates a long landing while a negative AT indi- 
cates a short landing. 

Figures 7, 8, and 9 show the descent of the DC -8 through the un- 
stable, neutral, and stable wind profiles of Figures 2, 3, and 4. The num- 
bering of the aircraft descent trajectories corresponds to the nvimbers on 
the wind profiles. The trajectories numbered 1 through 7 correspond to 
headwind profiles. The No. 8 trajectory corresponds to a zero wind (or 
constant wind) profile. The No. 9 and 10 trajectories correspond to tail- 
wind profiles. In particular, the No. 9 profile has the same shape as the 
No. 6 profile but differs in direction by 180 degrees. Similarly, the No. 10 
profile is the same as the No. 3 profile except for direction. The above is 
true for the unstable, stable, and neutral wind profiles. Figure 10 shows the 
DC-8 landing through the very stable wind profiles of Figure 11. The No. 1 
through 8 profiles are headwinds. No. 9 is a zero wind, and Profile No. 10 
is the same as No. 5 but is a tailwind. 
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Headwind Landings 


For headwind landings, the unstable wind profiles show the least 
scatter in touchdown point. As the stability increases, the aircraft be- 
comes more effected by the wind profile causing larger deviations in 
touchdown point. The very stable profiles produce very large deviations 
from the desired touchdown point, in one case over 3600 feet. In addition, 
the actual trajectory of the aircraft follows a steeper slope as the stability 
increases. The magnitude of the deviations in touchdown point for the 
stable and very stable conditions are significant and could lead to hazar- 
dous landing conditions. This is especially true since the simulation be- 
gan at an altitude of only 300 feet. In the very stable case, if the shear 
layer occurred above 300 feet, the effect on touchdown would certainly be 
greater than that shown in Figure 10. The same is not true for the other 
stability conditions. For the unstable, neutral, and stable profiles, the 
shear is close to the ground and little additional effect on touchdown point 
would result from beginning the simulation above 300 feet. Consequently, 
the very stable conditio n has the greatest potential for adversely affecting 
the landing of aircraft. 


Tailwind Landings 

For a tailwind, a somewhat larger deviation in touchdown occurs 
over the case of the same wind profile being encountered in the headwind 
direction. However, a landing accident involving a tailwind condition is 
probably not as likely to occur as one involving a headwind condition for two 
reasons. First of all, landings are made in the runway direction that has a 
headwind component whenever possible. As a result, tailwind landings, es- 
pecially with high wind magnitudes, are seldom required. Secondly, since 
a tailwind implies overshoot, the pilot has a slower descent rate and can 
more easily abort the landing. In addition, a light tailwind shear can pro- 
vide the somewhat desirable effect of a natural flare maneuver. 
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At this point, it is well to recall that the simulation model assumes 
no further control of the aircraft after trimming at the initial altitude of 300 
feet for a constant wind field; that is, no pilot feedback or automatic landing 
system is introduced. The deviation in touchdown point can therefore be con- 
sidered as a worst-case analysis in that any pilot or autopilot feedback during 
descent would, hopefully, result in smaller deviations in touchdown point. 

From this point of view, the Landing Simulation Program is intended to serve 
as a standard which will indicate the areas in which further research is required. 

Ground Effects 

In addition to wind shear causing an aircraft to depart from its desired 
2.7-degree glide slope, the ground effects on the aircraft aerodynamics near 
the surface will cause a small departure from the glide slope. To estimate the 
deviation in touchdown due to ground effects, simulations of the DC -8 and 
Boeing -747 landings were made with and without the inclusion of the ground 
effects terms for the wind profiles of Figures 2, 3, 4, and 11. Table 2 shows 
the results for the extreme stability conditions, unstable and very stable. For 
headwinds, the deviation in touchdown point is nearly identical with and without 
the ground effects included. For the zero wind case, the change in touchdown 
point is the largest but still less than 100 feet. For the tailwind cases in un- 
stable air, larger changes in touchdown point occur; however, this is somewhat 
misleading. In Case 6 of the Boeing -747, for example, the ground effects 
caused the aircraft to touch down at AT = 891 feet; whereas, without ground 
effects, the aircraft descended to an altitude of three feet at AT = 1320 feet 
downrange, then began to rise, finally touching down at AT = 3779 feet. With 
the pilot control, the actual difference between touchdown points would be 
considerably less. Consequently, since ground effects cause rather small 
deviations in touchdown and since accurate ground effect terms were not avail- 
able for all the aircraft used in this study, it was decided not to use the 
ground effects term in any further analyses of conventional aircraft. 
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TABLE 2 


DEVIATION FROM TOUCHDOWN POINT 
WITH AND WITHOUT GROUND EFFECTS 


Wind 

Profiles 

B 

-747 

DC 



With 

Ground 

Effects 

W ithout 
G round 
Effects 

With 

Ground 

Effects 

Without 
G round 
Effects 

Zero Wind 

- 90 

0 

- 91 

0 

No. 1 

- 566 

- 502 

- 553 

- 492 

No. 2 

- 852 

- 800 

- 834 

- 786 

No. 3 

- 821 

- 783 

- 809 

- 767 

^ No. 4 

+ 738 

+1049 

+ 770 

+1033 

rt No. 5 

- 890 

- 881 

- 888 

- 858 

« No. 6 

+ 891 

+3779 

+ 996 

+ 3715 

§ No. 7 

-1314 

-1317 

-1297 

-1281 

No. 8 

- 859 

- 858 

- 862 

- 835 

No. 9 

-2651 

-2652 

-2580 

-2579 

No. 1 

-3671 

-3677 

-3668 

-3669 

No. 2 

+ 1854 

+ 1833 

+2161 

+2167 

4 ) No. 3 

-2048 

-2065 

-2057 

-2056 

3 No. 4 

-2685 

-2650 

-2625 

-2625 

^ No. 5 

-2992 

-2965 

-2932 

-2933 

tK No. 6 

-2758 

-2762 

-2734 

-2735 

« No. 7 

-1641 

-1644 

-1619 

-1619 

^ No. 8 

- 715 

- 709 

- 675 

- 668 

No. 9 

- 333 

- 299 

- 305 

- 265 
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Wind Shear Effects on Touchdown 


The equation used to define the wind profile for unstable, neutral, 
and stable atmospheric conditions is defined in terms of the parameters 
u*, and L. The very stable profiles use as parameters the wind 

magnitude above the shear level and the altitude at which the shear occurs. 
Landing simulations of the DC -8 have been made to determine the influ- 
ence of each of these parameters on touchdown point for headwind landings. 
Figure 12 shows the results for the unstable case. For ease of interpre- 
tation, the X-axis does not use the frictional velocity, u*, but rather the 
wind speed at the initial height of 300 feet. From Figure 12, it is observed 
that the stability length, L, has little effect on touchdown point. The 
larger negative values of L, which imply increased stability, produce 
somewhat larger deviations in touchdown point. The surface roughness 
length, Z^, however, causes a larger variation in touchdown point. Over 
a very smooth surface such as mown grass (Z^ « 0. 001m) deviation in 
touchdown point, AT, greater than 1000 feet cannot occur even under a 
very strong wind field. If, however, the runway is surrounded by large 
buildings, with associated roughness length of one to three meters, AT 
increases by a factor of 2 to 3. 

The neutral wind profiles. Figure 13, for headwind landings can 
be considered as the limiting case of the unstable profiles as L - <x> . 

Thus, for a given Z^, AT will be larger for the neutral than for the un- 
stable cases. The effect of Z^ on AT is, however, somewhat less in 
neutral stability. That is, for a given wind speed, the difference between 
AT for say Z =3m and AT for Z =lm is smaller for neutral than for the 
unstable condition. As stability further increases to the positive side, 
becomes less important while L exerts greater influence on AT 
(Figure 14). Under highly stable conditions, with say L = 10, very large 
values of AT are produced. The terrain roughness has almost no influ- 
ence on these values of AT. It is under these highly stable conditions that 
hazardous landing situations are most likely to occur. 
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Deviation In Touchdown Point, 



Figure 12. Deviation in Touchdown for DC -8 in 
Unstable Wind Profiles: Headwind. 


u at 300 ft. 


Figure 13. Deviation in Touchdown for DC -8 in Neutral 
Wind Profiles: Headwind. 







The values of AT for very stable profiles are shown in Figure 15 
for headwind landings as a function of the altitude at which the shear, or 
interface layer occurs. The very st able profiles cause a larger deviation 
in touchdown p oint than the other wind profiles, especially when the shear 
for the very stable profile occurs at a high altitude. It is also the very 
stable profiles that cannot be predicted by an analytical exp.ression in- 
volving meteorological parameters. Thus, during very stable atmospheric 
conditions with Richardson number greater than 0. 2, careful attention should 
be payed to the natural environment. 


Comparison of Different Types of Conventional Aircraft 

Figures 16 through 19 compare the results from the DC -8 with that 
of other aircraft. The aircraft considered span a large range of size, 
weight, and landing speeds. Table 1 can be consulted for the exact 
aircraft weight and landing speeds used in the simulation. In general, 
the variation in touchdown point, AT, is not largely dependent on the type 
aircraft. Even the C-130 which is slower and lighter performs very simi- 
larly to the other aircraft. It is interesting to observe that the relative 
sensitivity of different aircraft to a given wind field depends upon the wind 
field. For example, a stable profile with = 0. 1 and L = 10 (Figure 18) 
shows the C-130 to produce the largest values of AT and the C-135A the 
smallest. However, by changing the value of L to L. = 300, the result is 
that the C-130 produces the smallest values of AT with the C-135A falling 
in the middle range. 

The variation in AT resulting from the simulated landing of different 
type aircraft is smaller than the variation in AT due to surface roughness 
and to stability. For practical considerations, the type of wind profile 
that is hazardous to one type of aircraft is hazardous to all types - at 
least within the range of aircraft discussed in this report. 
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Variation in AT due to Aircraft Weight 


In simulating the landing of an aircraft, the gross weight of the 
aircraft was chosen approximately midrange between its tolerance ex- 
tremes. To determine the influence of the gross weight of an aircraft 
on AT, landings were simulated for the Boeing-747 with gross weights of 
400,000 pounds and 550, 000 pounds. These values correspond approxi- 
mately to minimum and maximum landing weight limitations for this air- 
craft. Figure 20 compares the values of AT for the very stable profiles. 
The difference between AT for the two gross weights is small; generally 
less than 200 feet. For the other stability cases, the difference between 
the AT's are even smaller and hence not presented in this report. To 
summa r ize, the gross weight of an aircraft appears to have little effect 
on touchdown point, especially when compared to the effects produced by 
surface roughness and stability. 

Variation in AT due to Cg Locations 

The location of the center of gravity, Cg, of an aircraft depends 
upon the gross weight and weight distribution of the aircraft. In most of 
the previous simulations, the Cg location was chosen approximately midway 
between the tolerance extremes. To observe the influence of the Cg loca- 
tion on touchdown point, simulated landings of the Boeing -747 were made 
with three different Cg locations; 15% Mean Ae rodynamic Chord; 25% 

Mean Aerodynamic Chord; and 33% Mean Aerodynamic Chord. The 15% 
and 33% Mean Aerodynamic Chord are the extreme allowable tolerances. 
Table 3 shows the results for the unstable and very stable profiles. The 
difference between AT values for the three Cg locations is negligible in 
all cases. The differences in AT for neutral and stable profiles were 
found to be intermediate between the unstable and very stable extremes. 
Thus, the Cg location seems to have little influence on the deviation in 
touchdown point due to wind shears. 


36 




Davlation lii Touchdowni AT 


W X 400, 000 lbs 
T 


Shear at 


= 550,000 lbs 80m 
T 


= 400, 000 lbs 
T 

= 550, 000 lbs 

Shear at 
60m 


. 400, 0001b 

W = 550,000 lb 

Shear at 
40m 



W^= 400,000 lb. 

W = 550, 000 lb. 

Shear at 


400,000 lb 
‘ 550,000 lb 

Shear at 


u at 300 feet 


Figure 20. Deviation in Touchdown for B-747 in Very Stable 
Wind Profiles with Different Landing Weights 
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TABLE 3 


DEVIATION FROM TOUCHDOWN POINT FOR VARIOUS 
LOCATIONS OF THE CENTER OF GRAVITY: 

B-747, Weight=400, 000 pounds 



Cg= 1 5% Mean 

Cg=25% Mean 

Cg=33% Mean 

Wind 

Aerodynamic 

Aerodynamic 

Aerodynamic 

Profiles 

Chord 

Chord 

Chord 

No. 1 

- 463 

- 462 

- 463 

No. 2 

- 737 

- 736 

- 737 

No. 3 

- 731 

- 731 

- 731 

« No. 4 

+ 951 

+ 942 

+ 942 

5 No. 5 

- 856 

- 855 

- 849 

tS No. 6 

+ 2877 

+ 2985 

+ 3633 

5 No. 7 

-1323 

-1317 

-1299 

No. 8 

- 849 

- 847 

- 837 

No. 9 

-2887 

-2845 

-2781 

No. 1 

-4028 

1 

-3953 

-3854 

No. 2 

+ 1164 

+ 1126 

+ 1260 

« No. 3 

-2370 

-2299 

-2204 

S No. 4 

-2936 

-2884 

-2821 

" No. 5 

-3212 

-3177 

-3138 

“ No. 6 

-3031 

-2985 

-2932 

S No. 7 

-1854 

-1827 

-1797 

> No. 8 

- 824 

- 818 

- 813 

No. 9 

- 352 

- 352 

- 351 
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ANALYSIS OF AUGMENTOR-WING STOL AIRCRAFT 


The landing of the augmentor-wing STOL aircraft described in the 
Aircraft Landing Model Section was simulated for the same unstable, neutral, 
stable, and very stable wind profiles identified in the Analysis of Conventional 
Aircraft Section. The trim conditions were defined for a 7-degree glide 
slope with descent again beginning at 300 feet. Under these conditions, the 
resulting touchdown point for a constant (or zero) wind field would be (ignoring 
for a moment the ground effects on the aircraft) 2443 feet downrange. Figures 
21 through 24 show the landing of the STOL aircraft in the unstable, neutral, 
stable, and very stable profiles previously presented as Figures 2, 3, 4, and 
11, respectively. For several of the tailwind profiles, in particular those 
for which the wind velocity exceeded 17 ft/sec, it was impossible to trim the 
aircraft to follow a 7-degree glide slope by controlling only the thrust magni- 
tude and elevator. The solution of the initialization subroutine showed a nega- 
tive thrust would be required for the STOL aircraft to descend the 7-degree 
glide slope at a constant velocity relative to the air of 118ft/sec. Thus, the 
results from those runs where negative thrust was assumed have been dis- 
carded. For several wind fields, the shear causes the aircraft to depart 
from the glide slope but the natural response of the aircraft eventually brings 
it back toward the glide slope. Trajectories No. 1 and 10 of Figure 23 and 
Trajectory No. 10 of Figure 24 are good examples of this phenomenon. These 
oscillations arise from the Phugoid mode of the aircraft. In several of the 
STOL trajectories, the aircraft was further from the glide slope at some point 
prior to touchdown than it was at touchdown. This was not true of the conven- 
tional aircraft flights. In the latter case, the period of the Phugoid oscillation 
was much larger so that touchdown occurred before the first quarter-cycle 
of the flightpath oscillation. 

For most of the STOL trajectories in which the above phenomenon 
was present, AT was not greatly less than the maximiom deviation from the 
glide slope during the flight. Nonetheless, the AT's observed for the STOL 
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Figure 21, STOL Descent Trajectories Through Unstable Wind Profiles. 



Figure 22, STOL Descent Trajectories Through Neutral Wind Profiles. 












scent Trajectories Through Very Stable Wind Profiles. 



aircraft cannot be interpreted in quite the same manner as they were for 
the conventional aircraft. The STOL AT's only approximate the maximum 
deviation from the glide slope for a given wind shear. 

Ground Effects of STOL Aircraft 

The landing of the augmentor-wing STOL aircraft was simulated 
both with and without ground effects terms to determine the degree to which 
the aircraft ground effects influenced touchdown point. Table 4 shows the 
deviation in touchdown for the very stable and unstable wind profiles for 
the simulations with and without ground effects. For the zero wind field, 
the ground effects causes a touchdown 180 feet beyond the 7-degree glide 
slope touchdown point. This contrasts with the 90-feet-short touchdown 
point for the Boeing-747 and DC-8 aircraft. For the unstable profiles, the 
ground effects causes a deviation in touchdown point on the order of 150 
to 200 feet. For the very stable profiles, the deviations in touchdown 
are much smaller. In general, the ground effects provide the largest 
influence when the wind shear is light and the aircraft's descent follows 
a path near its desired glide slope. Since, even under these conditions, the 
influence of ground effects on touchdown is generally less than 200 feet, 
ground effects are not an important consideration in our study. 

Wind Shear Effects on Touchdown 

Figure 25 shows deviation in touchdown, AT, for selected stable 
and unstable wind profiles when the landing approach is in the headwind 
direction. For the stable wind profiles, AT shows the most variation with 
the parameter L. For a fixed value of L, the variation of AT with is 
small (not shown in Figure 25). For the unstable profiles, the variation 
in AT is greater with respect to than with respect to L. This is con- 
sistent with the previously observed variations in AT with Z^ and L for 
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TABLE 4 


DEVIATION FROM TOUCHDOWN POINT WITH 
AND WITHOUT GROUND EFFECTS 
FOR STOL AIRCRAFT 




Wind 

Profiles 

Augmentor-Wing STOL 

AT With 
Ground 
Effects 

AT Without 
G round 
Effects 

Difference 

Zero Wind 

+180 

0 

180 

No. 1 

+ 125 

- 91 

216 

No. 2 

+ 32 

-161 

193 

No. 3 

+ 57 

-145 

202 

iJNo. 4 

+349 

+190 

159 

(d No. 5 

- 8 

-179 

171 

«No. 6 

Negative T 

tirust Needed t 

> Trim 

DNo. 7 

-113 

-280 

167 

No. 8 

- 49 

-189 

140 

No. 9 

-599 

-713 

114 

No. 1 

-1299 

-1303 

4 

No. 2 

+ 157 

+ 137 

20 

No. 3 

- 683 

- 701 

18 

•2 No. 4 

wO 

- 975 

- 978 

3 

rtNo. 5 

-1141 

-1141 

0 

"^No. 6 

-1024 

-1025 

1 

^No. 7 

- 623 

- 639 

16 

>No. 8 

- 189 

- 290 

101 

No. 9 

+ 54 

- 123 

177 
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Figure 25. Deviation in Touchdown for STOL in 

Unstable and Stable Wind Profiles: Headwind. 



conventional -type aircraft. The variation of AT in neutral wind profiles, though 
not shown in Figure 25, lies intermediate between the results for the stable 
and unstable wind profiles. Figure 26 shows the variations in AT for the very 
stable wind profiles. The values of AT are approximately three times smaller 
for the STOL aircraft than for the conventional aircraft. In comparing the 
I STOL results to that of the conventional aircraft, it must be remembered that 
j the simulated landings for the STOL aircraft were down a desired glide slope 
of 7 degrees as compared to a 2. 7-degree glide slope for the conventional air- 
craift. The difference in glide slope angle allowed the STOL aircraft to descend 
i the glide slope in a shorter length of time which apparently accounted for the 
much smaller deviations in touchdown observed with the STOL. On the other 
hand, since a 2. 7 -degree slide slope is typical for a conventional aircraft, and 
7 degrees typical for a STOL, the comparison between the two is valid when 

! 

I considering the effect of wind shear on typical landing conditions at a given airfield. 

Landing of the augmentor-wing STOL aircraft in a tailwind presents 
a problem if the tailwind is greater than 17 ft/sec. With a 17 ft/sec tailwind, 

' zero thrust is required to maintain the 7-degree glide slope with a relative 
air velocity of 118 ft/sec. For a tailwind in excess of 17 ft/sec the seven- 
degree glide slope cannot be maintained without increasing the speed of the 
j aircraft as it descends. If it becomes necessary to land this STOL aircraft 
I in a large tailwind, the pilot should probably decrease his glide slope angle 
to two or three degree s if possible. 


i 
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Figure 26. Deviation, in Touchdown for STOL in 

Very Stable Wind Profiles: Headwind. 


SUMMARY AND CONCLUSIONS 


The three -degrees -of -freedom aircraft landing simulation study has 
determined the types of wind shear profiles that can produce potentially 
hazardous landing conditions. No pilot or auto-pilot feedback was introduced 
into the simulation. Deviations in touchdown point in excess of 3000 feet re- 
sulting from variation of the horizontal wind during the final 300 feet of 
descent have been observed under wind shear conditions that are not unrealistic. 
The influence of ground effects, center of gravity location, and gross weight 
of the aircraft on the deviation in touchdown point due to wind shears has also 
been investiaged. The specific conclusions resulting from this study are: 

a) Stable (0<Ri<0, 2) and very stable (Ri>0. 2) conditions are 
most likely to produce hazardous landing conditions. Deviations 
in touchdown of 2000 to 4000 feet have been observed for conven- 
tional aircraft. Neutral and unstable wind profiles seldom cause 
deviations in touchdown point in excess of 2000 feet for conven- 
tional aircraft and 600 feet for the augmentor-wind STOL. 

b) The deviation in touchdown point, AT, is more dependent upon 
the terrain roughness, Z^, than upon the stability length, L, 
under unstable and neutral wind conditions. Under stable wind 
conditions, the reverse is true. 

c) For very stable conditions, AT is most dependent upon the 
altitude at which the shear layer occurs. Very stable wind 
profiles are highly unpredictable and are not dependent upon 
the surface parameters (Z^, u*, L) that characterize the other 
stability conditions, 

d) For the aircraft considered here, the variation in touchdown 
due to the ground effects on the aircraft is small in comparison 

to the variation which can result from the wind shears investigated. 
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e) For a given wind profile, a tailwind direction produces a 
slightly larger deviation in touchdown than does a headwind 
direction. 

f) For the conventional airc raft, the size, type, and the landing 
speed of the aircraft has some influence on AT but this influence 
is considerably less than that due to surface roughness and 
stability length. The types of conventional aircraft studied were 
theC-130E, C-135A, C-141, DC-8, andB-747. 

g) The landing of the augmentor-wing STOL in a given wind field 
produced a much smaller value of AT than the landing of con- 
ventional aircraft in the same wind field. In particular, AT 
values for a STOL landing were 3 to 6 times smaller than the 
corresponding AT values for a conventional landing. The large 
glide slope angle for the STOL (7 degrees) allowed it to land 

in less time and was a major reason why AT was smaller for the 
STOL. Under all but the most extreme shear conditions, values 
of AT for the STOL aircraft did not exceed 1000 feet. 

h) The difference in touchdown points between the landing of a fully 
loaded aircraft and an empty aircraft was found to be small. The 
aircraft analyzed was the Boeing 747. 

i) The difference in touchdown point resulting from a shift in the 
Cg, within operational tolerance, was found to be negligible. The 
aircraft analyzed was the Boeing 747. In particular, this indicates 
that the Boeing 747 was well designed with respect to its response 
to wind shear. 

Aeronautical Safety Considerations 

This research program has provided results that have direct application 
to aeronautical safety at airports. It has been shown that the most critical 
landing conditions are likely to occur under the stable and very stable 
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atmospheric conditions. Unfortunately, the wind profile under very stable 
conditions cannot be determined from a single wind measurement knowing the 
surface roughness and stability parameters. Very stable conditions often occur 
during the night under strong temperature inversions. Under these conditions, 
height and magnitude of the shear can only be determined by empirical measure- 
ments. Thus, an important consideration in improving airport safety is to 
provide the capability of measuring a vertical profile of shear. A first-order 
shear approximation could be derived by a simple two-point shear, calculated 
from an airport surface wind measurement and an onboard aircraft wind speed 
measurement. A more refined profile could perhaps be derived by some of the 
remote sensing techniques presently under development. 

A second safety consideration relevant to air traffic control is the obser- 
vation that a severe shear condition effects all aircraft (at least those within 
the range of size and type considered in this report) to approximately the same 
degree. Thus, if one aircraft reports landing difficulties due to wind shear, 
it should be assumed that all -size aircraft will experience similar landing 
problems. 

A third consideration related to airport safety concerns the homogeneous 
terrain that immediately surrounds different airports. Since the surface 
roughness parameter, Z^, shows little influence on touchdown points under stable 
and very stable conditions, the roughness of the terrain surrounding the air- 
port is not an important consideration for determining which airports experience 
the critical shear profiles. This statement needs some qualification, however, 
since the wind profiles defined in this study assumed the terrain to be homo- 
geneous of a given roughness length. At many airports, as the aircraft descends 
the glide slope, the roughness of the terrain that regulates the profile changes 
so that more shear could possibly be introduced into the wind profile. A study 
should be directed toward the change in terrain problem before a firm con- 
clusion delineating the terrain effects of airport surroundings can be made. 
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A fourth safety consideration is related to pilot training and automatic 
landing systems. A wide variety of wind profiles representing the various terrain 
roughness and stability conditions should be used as input to flight simulators 
where pilot response or pilot training is required. The same variety of wind 
profiles should be used for the evaluation of an automatic landing system. 

A final safety consideration has application to those involved in STOL 
airport design. The deviations in touchdown observed for the STOL aircraft 
should provide preliminary guidelines for the runway length safety factor 
needed to allow for touchdown dispersions due to wind shear. 

Suggestions for Additional Research 

This study has provided some basic results concerning the amount of 
dispersion in touchdown that is likely to occur for an aircraft landing under a 
specific set of assumptions. These assumptions limit the reality of the simu- 
lation model to some extent since the actual situation is different when a pilot- 
controlled aircraft descends the glide slope. Many factors which were not in 
eluded in this study may be important and could be considered in a more defined 
six-degree-of-freedom aircraft simulation model. A six-degree -of -freedom 
model could extend the analysis to consider cross wind landings and landings 
in wind fields where the direction of the wind changes significantly over the 
final 300 feet of descent. The effect of turbulence could also be studied by 
introducing three-dimensional turbulence structure into the mean wind profile. 
Further work is also needed in defining the wind profile along the glide slope 
when a change in terrain roughness occurs. The effect of the roughness change 
from sea to land terrain should also be studied. A final recommendation is to 
improve the reality of the simulation by introducing pilot feedback or an auto- 
matic landing system into the simulation model. Only by introducing control 
can a completely realistic simulation be anticipated. Such a research program 
as defined above would produce more definitive results under a much broader 
range of conditions. 
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APPENDIX A 


The equations for the lift, drag, and moment coefficients (C_ ,C ,C ) 

L D m 

and the aerodynamic coefficient data for each aircraft considered in this study 
are presented below. The ground effects terms for four of the aircraft (DC -8, 
C-135A, Boeing 747, and the STOL) are discussed. Ground effects data were 
not available for the C-130E and the C-141 aircraft. The aerodynamic data 
presented corresponds to a given flap setting and, where applicable, a given 
horizontal tail setting. The data includes the effects of the landing gear. 


DC -8 

The expressions for the aerodynamic coefficients of the DC -8 
aircraft are: 


_ c a _ c a 

Ct=C_ +C- a’+C-Te + ->ir ACv 


GE 


^D = CDq + Cd<i ® ' + ^^Dge 


C_ = C_ o' + C 


m 


o ntio 


c q ca ^ 

E ^^a ™GE 


where the usual notation is used for the various stability derivatives. The 
last term in each of the above equations is the ground effects term. The 
values of the stability derivatives are given in Table A. 1 for a flap setting 
of 50 degrees. The ground effects terms for the DC -8 are: 


AC 

AC 

AC 


lge 

I>GE 

“GE 


= 0.063 

= (-0.02 - 0.332q') Sqj. 


= - 0.066 


A-1 


TABLE A. 1 


DC -8 AERODYNAMIC DATA 


C, 

0. 90 




5. 30/ rad 

^L6e 

0. 0053/deg 

CLq 

7. 68/rad 

^La 

0. 0 

^Do 

0. 140 


0. 501/rad 

^Da2 

1. 818/rad^ 

^mo 

-0. 01 

^ma 

-1. 062/rad 

Crngj, 

-0. 0161/deg 

Cniq 

-12, 30/rad 

Cmi 

-4. 01 /rad 


A-2 


where 


'gE = 0.972e 


-h/17 


h is the wheel height, and C is the free stream value of C 


oo 


C-135A 

The expressions for C, and C for the C-135A are the same as 
^ L m 

those for the DC -8. The expression for the drag coefficient is 

Lj 

Table A. 2 shows the aerodynamic data for a flap setting of 30 
degrees and a horizontal tail deflection of -4 degrees. The ground effects 
terms are: 

AClqe = 039 + 0. 2292a •) 

^^Dge = ^4o ^ ®GE 

(0. 0228 - 0. 1408 Ct + 0. 054 0^ ) 

niGE L L ' GE 

CD CD 

where 

^ -h/29 

e = e 


C-141 


The expressions for the lift, drag, and pitching moment coefficients 
for the C-141 are the same as those for the C-135A. Table A. 3 presents 
the aerodynamic data for a flap setting of 45 degrees and a horizontal tail 
setting of -6 degrees. No data were available for C-141 ground effects 
terms. 
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TABLE A. 2 


C-135A AERODYNAMIC DATA 


Cl 

0.612 

O 

Cl 

a 

4 .01 /rad 

<=L6e 

0. 00376/deg 


0. 0 


0. 0 


0. 0685 

C i“v 2 

^ct 

0.0473 


0. 0922 


-0. 765/rad 

^msE 

-0. 0108/deg 

Cmq 

-14. 182/rad 

Cmi 

-5.787/rad 
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TABLE A. 3 


C-141 AERODYNAMIC DATA 



1. 309 

•"La 

5. 441 rad 


0. 00443/deg 

Cl, 

0. 0 

Cl 

0. 0 

^Do 

0. 0835 

CdcI 

0. 0388 

Cm© 

0.391 


-1. 351/rad 

^“6e 

-0. 0149/deg 

Cm„ 

-15. 75/ rad 

Crr^- 

-5. 17/rad 
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C-130E 


The C-130E aircraft is powered by four propjet engines. For this 
type of aircraft several of the stability derivatives depend upon the thrust 
coefficient defined as 

'^C = 

where is the thrust per engine, q is the dynamic pressure, and d 

is the propeller diameter. For the C-130E, d is equal to 13. 5 feet. The 
expressions for Cj_^, and for the C-130E are the same as for the 

C-135A except that ^mQ» upon T^. This 

dependence is shown in Table A. 4. The data in Table A. 4 is for a flap 
setting of 18 degrees. No ground effects data were available for this aircraft. 


Boeing 747 

The expressions for C_ , C-r\, and C used in this study for the 

Ij D m 

Boeing 747 are the same as those of the DC -8 with the exception that an 

additional term has been added to C to account for different center of 

m 

gravity locations. This term is: 

Cl ( Cg - 0.25) 

where Cg is the location of the center of gravity in terms of the mean 
aerodynamic chord, c . For a center of gravity location of 25% c, the 
above term is zero. The aerodynamic data presented in Table A. 5 is for 
a flap setting of 30 degrees and a horizontal tail setting of -4 degrees. 

It was assumed that the inboard and outboard elevator deflections are the 
same. The expressions for the ground effects terms are: 


TABLE A. 4 


C-130E AERODYNAMIC DATA 



0.379 T^ + 0. 702 


(3.62 T^ + 6. 70)/ rad 


0. 0092/deg 


6. 59/ rad 

CLa 

2. 52/ rad 

^Do 

0. 0638 

^Dcl 

0. 0305 

Cr^ 

hIq 

-0. 224 Tc + 0- 338 


(2.75 T^ - 1.785)/ rad 


-0. 0285/deg 

niq 

-20. 06/ rad 

c„« 

“Vi 

-7. 91 /rad 
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TABLE A 5 


BOEING 747 AERODYNAMIC DATA 


0 

U 

0. 960 


5. 735/ rad 


0. 0062/deg 

Cl 

5. 68/ rad 
6. 76/ rad 
4. 88/ rad 


-6. 70/ rad 

^Do 

0. 1381 


0. 5498/ rad 


2. 190/rad^ 

Cmo 

0. 094 

Cm 

-1. 536/rad 


-0. 0246/deg 


-21. 50/ rad 


-3. 40/rad 
-3. 81 / rad 
-3. 09/rad 


for c. g. 
for c. g. 
for c. g. 


for c. g. 
for c. g. 
for c. g. 


= 25% 
= 15 % 
= 33% 


25% c 
1 5% c 
33% c 


A-8 


lul u I o 



= K® (0. 240) cos rs. 036(a '-0. 00526)] 

J-*GE GE 

AC^ = K;^^ (2. 308 a - 0. 9796a 0. 1769o ' - 0. 0384) 

°GE 

*^GE - 0-115) 


where 


= 1. 7034 X 10"^ - 1. 0736 x lO""* - 1. 4813 x lO"^ h + 1. 0 

GE 

K® = 3. 7906 X lO'^ - 4. 937 x lO'^ + 2. 807 x lO"^ h + 1. 0 
GE 


The ground effects terms are included only during the last 82. 5 feet of flight. 


STOL 


The equations for the aerodynamic coefficients for the augmentor 
wing STOL aircraft are somewhat different than for the previously discussed 
aircraft. Several of the stability derivatives are functions of the thrust co- 
efficient Cj. In the present study, a value for of 0. 75 was chosen. The 
aerodynamic coefficients are expressed as: 


C_ = Ct +C, +C,„ a'+Cr- 6 f+^Ct C. 

^ ^WB l^Ho E 2V^ 


Cd = CDo + Cu^a • + 

^LWB +-^ '‘"o ’> 

c c 


+ sino ' - ^ cosa ') + (C^jj + a ’ +C^c Se> 

c c o a jE 

sina' + AC 


oige 
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The subscripts H and WB refer to the contributions from the horizontal 

tail and wing-body, respectively. C, is the basic lift coefficient for 

^WB 

the wing-body. The parameters Z-y^, and are distances which 
relate the aircraft center of gravity to the wing -body reference point 
and the 25% mean aerodynamic chord of the horizontal tail. (See Figure A. 1). 
The values of these parameters are: 


fyy. = 3. 95 ft. 
Zyy = 0. 083 ft. 
f jj = 75. 39 ft. 
Zj^ = -24. 75 ft. 


The values for the various stability derivatives are shown in Table A. 6. 
The flap setting is 70 degrees and the auxiliary flap is set at 6 degrees. 


The description of the ground effects for the STCL. aircraft is much 
more complex than for the conventional aircraft investigated. The ground 
effects come into play at approximately 200 feet altitude. Prior to this 
point, is given by 

^Lwb^"" ^LWBq ^ ^LwBa ' 


The subscript oo denotes the free stream value. The value of C 
modified by ground effects according to the equation 


l-wb 


is 


oo 


"WB 


GE 




whe re 


q/q. 


C30 


c 


1+ LwEbo/___l \ 


c • Ar 
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eters for the 
Aire raft. 


TABLE A. 6 


AUGMENTOR-WING STOL AERODYNAMIC DATA 



4.274 


4. 498/ rad 


-0. 124 


0. 743/rad 


0. 00936/deg 


0. 0 

a 

0. 0 

^Do 

0.449 


2. 12/ rad 

^nrio 

-1.245 

Cm, 

0. 372/ rad 

Cmq 

-0. 0265/rad 

Cmi 

0. 0 

Clc. 

2.20 

J 


CDc. 

0. 1709a' - 0 


-0. 780 

Cin 

0. 384/ rad 


a 
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ZttAjj^ [1 + 16 ( 


c • A, 




and A— is the wing aspect ratio which has a value here of 6. 5. The above 
equation does not apply if C- exceeds its maxim\im allowable value 
which for the present case is 6. 74. The drag and pitching moment ground 
effects terms are: 




CO 


+ (1 . q/q ) C Cj 
Qo Un • J 

'J 


AC 


An 


mCE 


= a -) (q/q -1) + C_ (-p- ) C 


CO 






+ ( 1- q/q ) 


Values for Cl^,. Cd^- . C_^(- 
for a Cj value of 0. 75, ^ 


, and C’ are presented in Table A. 6 


m„ 
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